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■ Abstract. We study the viability of the lightest neutralino as a dark matter candidate in the Next- 

to-Minimal Supersymmetric Standard Model. Taking into account accelerator constraints as well as 
bounds on low-energy observables, and imposing consistency with present bounds on the neutralino 
relic density, we address the prospects for the direct detection of neutralino dark matter. We find 
regions of the allowed parameter space where the neutralino detection cross section is within the 
^j. . reach of dark matter detectors, essentially owing to the presence of very light singlet-like Higgses, 

' and to either singlino dominated or very light neutralinos. 
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Ph| The Next-to-Minimal Supersymmetric Standard Model (NMSSM) is a well- 

motivated extension of the Minimal Supersymmetric Standard Model (MSSM) by 
a singlet superfield S. In addition to providing an elegant solution to the so-called /! 
problem of the MSSM, and rendering the Higgs "little fine tuning problem" of the 
> ' MSSM less severe, the presence of additional fields (CP-even, CP-odd neutral Higgs 

bosons, and a fifth neutralino) leads to a richer and more complex phenomenology. In 
^ ! the NMSSM, one has possibility of dark matter scenarios that can be very different from 

those encountered in the MSSM, both regarding the relic density as well as prospects 
for direct detection. In particular, the exchange of very light Higgses can lead to large 
direct detection cross sections, within the reach of the present generation of dark matter 
detectors 01,0]. 

The NMSSM superpotential contains new couplings, involving the Higgs doublets 
•£j . and the singlet, 

X: 

Wnmssm = Bij (y u H{Q 1 u + Y d H[QJ d + Y e H\V i) - £,jXSH[H^ + ^kS 3 . (1) 

New terms associated with additional soft supersymmetry (SUSY) breaking trilinear 
couplings Ax and A K will also appear in the Lagrangian. After the spontaneous breaking 
of electroweak (EW) symmetry, the neutral Higgs scalars develop vacuum expectation 
values (VEVs), (H®) = v\, (Z/^) = v 2 an d (S) = s. This leads to the dynamical generation 
of an effective interaction 11H1H2, with /! = Xs. In the NMSSM spectrum, we now 
have three CP-even and two CP-odd Higgs states. In particular, the lightest Higgs 
scalar can be written as h® = Si\H® + SiiH® + S13S, where S is the unitary matrix that 
diagonalises the 3x3 scalar Higgs mass matrix. In the neutralino sector, the singlino 
mixes with the bino, wino and Higgsinos. The lightest state can be now expressed 
as Xi = N n B° + Ni 2 W$ + N U H^ + N U H% + N l5 S, where N diagonalises the 5 x 5 
neutralino mass matrix. 
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The low-energy NMSSM parameter space can be described in terms of the 
A, fc, tan/3, /i, Ax, A K degrees of freedom, as well as the soft SUSY-breaking terms, 
namely gaugino masses, M\ 2,3. A thorough analysis of the low-energy NMSSM 
phenomenology (minimisation of the potential, computation of spectrum and compati- 
bility with LEP/Tevatron bounds) can be obtained using the nmhdecay 2.0 code (Hi. 
Additionally, we have also included in our analysis [2] a more precise computation 
of the b — > sy decay in the NMSSM Q4£| , taking into account next-to-leading order 
contributions, and imposing consistency at the 2a level with the experimental central 
value [H,BR ex P(Z?^y) = (3.55 ±0.27) x 10 4 . Likewise, we have also incorporated 
the constraints coming from the contribution of a light pseudoscalar a in the NMSSM 
to the rare B- and K- meson decays [|4J]. Finally, in our analysis we have also included 
the constraints coming from the SUSY contributions to the muon anomalous magnetic 
moment, = (g^ — 2). At present, the observed excess in a^ xp [6] constrains a possible 
SUSY contribution to be [0] aJ USY = (27.6 ± 8) x lO" 10 . 

As thoroughly discussed in (TLB], an extensive part of the low-energy NMSSM pa- 
rameter space is directly excluded on theoretical grounds, namely the occurrence of 
tachyons. Moreover, false minima, and Landau poles for the couplings in Wnmssm 
also render unviable further areas. In regions surviving the latter constraints and the 
LEP/Tevatron bounds, and that present the most appealing prospects regarding dark mat- 
ter direct detection, an important role is played by the b — > sy decay, which can in prin- 
ciple exclude important regions of the parameter space (due to very large contributions 
to BR(b — ► sy) from charged Higgs exchange). Concerning the SUSY contributions to 
dp, in the dark-matter interesting regions, these tend to be in general quite small. A suf- 
ficiently large a^ USY can nevertheless be obtained when slepton (and gaugino) masses 
are decreased, in association with large values of the slepton trilinear couplings. 

In order to be a good dark matter candidate, the lightest NMSSM neutralino must also 
comply with the increasingly stringent bounds on its relic density, 

0.1 < Qk 1 < 0.3 (astrophysical), 0.095 < Qk 1 < 0.112 (WMAP), (2) 

respectively arising from astrophysical constraints J[8j], and from taking into account the 
recent three years data from the WMAP satellite I9J]. 

Compared to the MSSM, there are several alterations leading to Q.^oh 2 : first, the 
possibility of a singlino-like lightest supersymmetric particle (LSP), associated with 
new couplings in the interaction Lagrangian, may favour the coupling of WIMPs to a 
singlet-like Higgs, whose mass can be substantially lighter than in the MSSM, given 
the more relaxed experimental constraints. Secondly, in the NMSSM we have new 
open channels for neutralino annihilation, e.g. s-channel resonances, new channels with 
annihilation into Zh®, h®h®, h®a® and a^a® (due to light h® and a® states), providing 
important contributions to the annihilation and co-annihilation cross-sections ft 1 011 - In 
our analysis |2D, the results for the neutralino relic density, obtained from an nmhdecay 
link to MicrOMEGAS 03] . 

In the regions of the parameter space likely to have large neutralino detection cross 
sections |Q]], we have found that, in general, the correct relic density can only be obtained 
when either the singlino composition of the neutralino is large enough or when the 
annihilation channels into Z, W, or h® are kinematically forbidden. Interestingly, some 



allowed areas were very close to the tachyonic border, which as we verified, can give 
rise to very large direct detection cross sections. 

As pointed out in [1], the existence of a fifth neutralino state, together with the pres- 
ence of new terms in the Higgs-neutralino-neutralino interaction (which are proportional 
to A and k), trigger new contributions to the spin-independent part of the neutralino- 
nucleon cross section, o„o_ . On the one hand, although the term associated with the 

K\ P 

^-channel squark exchange is formally identical to the MSSM case, it can be signifi- 
cantly reduced if the lightest neutralino has a major singlino composition. On the other 
hand, and more importantly, the dominant contribution to a^o _ p , associated to the ex- 
change of CP-even Higgs bosons on the ^-channel, can be largely enhanced when these 
are very light. Consequently, large detection cross sections can be obtained, even within 
the reach of the present generation of dark matter detectors. 

As illustrative examples of our analysis, let us consider two cases, which represent the 
most relevant features of NMSSM dark matter scenarios. Let us begin by considering 
Mi = 160 GeV, A x = 400 Ge V, A K = -200 GeV, and \i = 130 GeV, with tan j8 = 5, lead- 
ing to a parameter space consistent with bounds on fl^ USY and BR(£> — > s y). In regions 
of the parameter space where the neutralino is relatively heavy and has a mixed bino- 
Higgsino composition, the relic density is too small. As the neutralino mass decreases, 
some annihilation channels become kinematically forbidden, such as annihilation into 
a pair of Z or W bosons when m^o < Mz or m^o < Mw, respectively. Below these, the 
resulting relic density can be large enough to fulfil the WMAP constraint, which occurs 
for distinct regions in parameter space B20. One of the allowed regions is close to the 
tachyonic area and exhibits very light singlet-like Higgses, potentially leading to large 
detection cross sections. This is indeed the case, as evidenced on the left-hand side of 
Fig. CD where the theoretical predictions for c^o_ p are plotted versus the lightest neu- 
tralino mass. The resulting G%o_ p spans several orders of magnitude, but, remarkably, 
areas with o^o > 10~ 7 pb are found. These correspond to the above mentioned regions 
of the parameter space with very light singlet-like Higgses (25 GeV < m h o < 50 GeV 

with S 2 3 > 0.99). The neutralino is a mixed singlino-Higgsino state (N 2 5 ~ 0.35) with 
mass around 75 GeV. The sensitivities of present and projected dark matter experiments 
are also depicted for comparison. On the right-hand side of Fig. Q] we show the resulting 
<Xy0 when the neutralino composition is changed, namely when the Higgsino com- 

ponent is enhanced. In particular, we consider the choice Mi = 330 GeV, tan/3 = 5, 
Ax = 570 GeV, A K = —60 GeV, with /i = 160 GeV. Such neutralinos annihilate more 
efficiently, thus leading to a reduced h 2 , so that the astrophysical constraint becomes 

more stringent. On the right-hand side of Fig. [H the various resonances appear as funnels 
in the predicted o^O-p f° r the regions with the correct £l^o h 2 at the corresponding values 
of the neutralino mass (m^o ~ Mz/2 and m^o ~ m h o/2). Below the resonance with the 
Z boson, light neutralinos are obtained m«o < Mz/2 with a large singlino composition 
which have the correct relic abundance. The lightest Higgs is also singlet-like and very 
light, leading to a very large detection cross section, o^o > 10~ 6 pb. 




FIGURE 1. Scatter plot of the scalar neutralino-nucleon cross section as a function of the lightest 
neutralino mass. On the left, M\ = 160 GeV, tan/3 = 5, Ax = 400 GeV, A K = -200 GeV, and jj, = 130 
GeV. All the points represented are in agreement with LEP/Tevatron, a^ USY , and BR(b — » sy) bounds. 
Dark grey dots represent points which, in addition, fulfil 0. 1 < Q. 9 o h 2 < 0.3, whereas black dots are those 
in agreement with the WMAP constraint. The sensitivities of present and projected experiments are also 
depicted, with solid and dashed lines, respectively. On the right, a different example with Mi = 330 GeV, 
tan/3 = 5, A x = 570 GeV, A K = -60 GeV, with ji = 160 GeV, a case where the resulting a^ USY is outside 
the experimental 2a region. 
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